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Abstract 

Ground -based spectrometers have been developed to measure the concentration, 
velocity, and temperature of ozone in the mesosphere and lower thermosphere (MLT) 
using low-cost satellite television electronics to observe the 11.072 GHz spectral line 
of ozone. A two-channel spectrometer has been engineered to yield various perfor¬ 
mance improvements, including a doubling of the signal-to-noise ratio, improved data 
processing efficiency, and lower power consumption at 15 W. Following 2009 and 2012 
observations of the seasonal and diurnal variations in ozone concentration near the 
mesopause, the ozone line was observed at an altitude near 95 km and latitude of 38 
degrees north using three single-channel spectrometers located at the MIT Haystack 
Observatory (Westford, MA), Chelmsford High School (Chelmsford, MA), and Union 
College (Schenectady, NY) pointed south at 8 degrees. Observations from 2009 through 
2014 are used to derive the nightly-averaged seasonal variation in meridional velocity, 
as well as the seasonally-averaged variation with local solar time. The results indicate 
a seasonal trend in which the winds at 95 km come from the north at about 10 ms -1 in 
the summer of the northern hemisphere, and from the south at about 10 ms -1 in the 
winter. Nighttime data from -5 to +5 hours local solar time show a gradual transition 
of the meridional wind velocity from about -20 ms -1 to +20 ms -1 . These two trends 
correlate with nighttime wind measurements from the Millstone Hill High-Resolution 
Fabry-Perot Interferometer (FPI) in Westford, MA, which uses the 557.7 nm green line 
nightglow from atomic oxygen centered at 95 km. The results have also been compared 
with average meridional winds measured with meteor radar. 
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1 Introduction 


1.1 Past and New Work 

The Mesospheric Ozone System for Atmospheric Investigations in the Classroom (MO¬ 
SAIC) is an array of single-channel radio telescopes that observe the 11.072 4545 GHz 
spectral line of ozone in the middle atmosphere. System hardware is based on the 
Very Small Radio Telescope (VSRT), which was developed at the Massachusetts Insti¬ 
tute of Technology (MIT) Haystack Observatory to perform radio observations of the 
Sun and demonstrate physical principles related to interferometry and electromagnetic 
radiation. Although the VSRT and single-channel MOSAIC systems were originally 
intended to be educational outreach tools, they have since been deployed to conduct 
original scientific research measuring the diurnal, seasonal, and latitudinal variations in 
ozone concentration [1] [2] . In recent years, it was realized that a greater resolution of 
the spectral line was needed in order to produce more reliable data about the ozone’s 
velocity and temperature. A six-channel unit was engineered in order to improve the 
signal-to-noise ratio by V6, and has been in operation for less than one year [3]. This 
technical report encompasses the engineering of an improved two-channel system [4]. 
It also presents a detailed study of ozone velocity and preliminary observations of ozone 
temperature as obtained from a network of single-channel spectrometers observing the 
spectral line at 95 km. Additionally, we compare our results with the 557.7 nm oxy¬ 
gen spectral line, as well as average winds measured with meteor radar and various 
theoretical and empirical models of the middle atmosphere. 

1.2 Ozone Dynamics in the Middle Atmosphere 

The ozone molecule consists of three oxygen atoms (O 3 ) and is a gas that thrives in 
various regions of the Earth’s atmosphere. The vast majority of ozone resides in the 
stratosphere, which is 10-50 km from the Earth’s surface, and forms the “ozone layer” 
that shields the Earth from solar ultraviolet radiation [5]. In much smaller quantities, 
ozone exists in the mesosphere, which is 50-80 km from the surface, as well as the 
mesopause, the region between 80-85 km and 100-105 km [ 6 ] that acts a minimum 
temperature boundary between the mesosphere and the thermosphere. Although the 
thermosphere itself is the region above 85 km [7], we shall refer to our target altitude of 
95 km as being “near the mesopause.” As a whole, the region of interest is denoted the 
“mesosphere and lower thermosphere” (MLT) [ 8 ] [9]. In this region, new ozone molecules 
are created through only one atmospheric chemical process: 

0 + 0 2 + M^0 3 + M, (1) 

which involves a collision between atomic oxygen (O), diatomic oxygen ( 0 2 ), and a third 
element M, which usually represents diatomic nitrogen (N 2 ) or diatomic oxygen. There 
are two processes that destroy ozone in the mesopause, the first being a collision with 
a single photon that has energy E = fno , where h = 1.055 x 10 ” 34 J • s is the reduced 
Planck constant, and ui is angular frequency in the ultraviolet band, 750 THz < ui/2n < 
3000 THz. This event splits ozone into atomic and diatomic oxygen: 

0 3 + hu;^ 0 + 0 2 . ( 2 ) 
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Figure 1: Daytime ultraviolet annihilation of ozone. 


The second destructive process is a collision with atomic hydrogen, which produces the 
hydroxyl free radical (OH) and diatomic oxygen: 

0 3 + H+> 0H + 0 2 . (3) 

During the daytime, almost all of the ozone in the mesopause is destroyed by ultraviolet 
photon collisions. This results in a significant decrease in ozone concentration, followed 
by a nighttime increase [1] [2] [10] [11]. This is shown in Figure 1, where the ozone volume 
mixing ratio (in parts per million) drops to zero between +6 and +18 hours local solar 
time. The data were taken from the ozone spectrometer located at MIT Haystack 
Observatory (42.618° N, 71.498° W) and averaged for the first 204 days of 2015. 

1.3 Atmospheric Winds and Tides 

As microwaves propagate at the speed of light, c, the Doppler shift of the spectral line 
from one frequency /o to another, 

f= ( 1 + l) /o ’ (4) 

directly produces the ozone velocity, v, which can reveal information about winds, tides, 
and other forms of wave propagation in the atmosphere. For example, gravity waves are 
the result of a gravitationally-driven buoyant restoring force on circulating tropospheric 
air masses. They dissipate in the stratosphere and mesosphere and heavily influence 
global air flow [9] [12]. Planetary Rossby waves are generated from temperature and 
pressure fields that arise from the Coriolis effect attributed to the Earth’s rotation 
[12]. However, in the MLT, atmospheric tides are most heavily dominant instead of 
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Figure 2: Two-channel spectrometer block diagram. 


gravity waves and planetary waves. The solar ultraviolet heating of stratospheric and 
mesospheric ozone makes atmospheric tides a primary factor in the diurnal trend in 
ozone velocity. Further importance arises from tidal momentum transfer dynamics, 
which have been known to be responsible for the coupling between the upper and lower 
atmosphere [8] [9]. Thus, studying the seasonal and temporal trends in atmospheric tides 
obtained from wind velocity data are critical to understanding the physics of global air 
circulation in the middle atmosphere. 


2 The Two—Channel Ozone Spectrometer 

2.1 Electrical Hardware and Software 

The upgraded two-channel MOSAIC hardware (Figure 2) first consists of a Winegard 
DS-4040 parabolic reflector with 45.72 cm diameter. It is tuned to the 11.072 4545 GHz 
ozone spectral line by means of a frequency calibration pulse generator, which uses a 
Connor Winfield OX200-SC 10 MHz oven crystal oscillator to generate a harmonic at 
11.070 GHz. It should be noted that the ozone line is in the X band of the microwave 
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Figure 3: Spectrum of the frequency calibration signal. 


spectrum [13] but can be resolved with commercially available K u band satellite tele¬ 
vision electronics, whose waveguides can accept frequencies at the far edge of the X 
band. 1 Specifically, we use a Star-Corn SR-3602 Mini low-noise block downconverter 
feed (LNBF) that covers 10.70-11.70 GHz on the low-band and 11.70-12.75 GHz on the 
high-band [14]. The LNBF has a 9.75 GHz local oscillator whose signal is heterodyned 
with the ozone line to downconvert to 1.32 GHz in the L band. The spectrum of the 
frequency calibration signal (Figure 3) is observed by tuning to a line frequency of 1.32 
GHz (corresponding to 11.070 GHz) and selecting a command to leave the frequency 
calibration on [4] . 

A 12 V wall-wart is used to select for the horizontal microwave polarization (H) on 
Channel 1, and an 18 V wall-wart is used to select the vertical microwave polarization 
(V) on Channel 2. Next, in order to minimize out-of-band interference, the signal from 
each channel is separately fed into a bias tee, followed by a Very High Frequency (VHF) 
high-pass filter (HPF) at 740 MHz, and finally an 8 dB signal attenuator. Each signal 
then enters an R820T RTL2832 USB TV dongle that has a local oscillator tuned to 
the range around 1.42 GHz in order to receive the ozone line. The ozone signal is 
heterodyned with the local oscillator and then sent to the computer for processing. 

Our choice of machine is Intel’s Next Unit of Computing (NUC), which is small form- 
factor device, similar to the Raspberry Pi and the BeagleBoard BeagleBone, but more 
powerful. It runs the Ubuntu Linux 14.04 LTS “Trusty Tahr” operating system. The 

1. This is especially true in Europe, where television signals are commonly broadcast at upper X band 
frequencies around 11 GHz. Thus, an ozone spectrometer in Europe would pick up television interference, 
while a spectrometer in the United States would be quiet since U.S. television signals are encoded in the 
K u band above 12 GHz. Because of this distinction, our spectrometers use LNBFs designed for European 
satellite television dishes. 
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source code in C performs a dual-core threaded Fast Fourier Transform (FFT) on the 
signal to extract the ozone spectrum, and sends the resulting data via Internet to the 
MOSAIC GUI at www.haystack.mit.edu/ozone/. For the physical layout of the 
two-channel spectrometer hardware, see Figure 7. 

The optimum digital signal levels after 8 dB attenuation are about mid-scale in the 
ADC range of -127 to +128. Frequency switching is used, in which the ozone line is 
alternately placed at 0.9375 and 1.5625 MHz in the baseband 2 output. While the line 
is at 0.9375 MHz, the spectrum at 1.5625 MHz is used as a reference, and vice versa. 
In addition, the reference spectrum is smoothed using a polynomial. In this case, the 
noise in the line spectrum is reduced so that the normalized theoretical noise is given 
by the inverse of the square root of the resolution bandwidth times the integration 
time as would be obtained by a total power spectrometer. When the spectra of both 
polarizations are combined, the theoretical noise level for 80 K system noise and 2.44 
kHz resolution is 19 mK in one hour. Software timing tests show that for 42 seconds of 
data (1.05 x 10 8 complex samples in each channel), the FFT and other calculations add 
5 seconds, and the calibration another 10 seconds. Most of the time for calibration is 
the result of using stty to turn on and off the RS232 RTS line. In principle, this added 
time could be reduced by using another device, such as a USB power switch dongle. 

2.2 Performance Improvements 

The combination of using both (H) and (V) polarizations and smoothing the reference 
spectra results in an improvement in the signal-to-noise ratio (SNR) by a factor of 2. In 
order to further improve the SNR, we selected for the best LNBF by characterizing the 
noise figures of five different LNBF models. A thorough description of our procedure, 
and complete results, can be found in Appendix §9. In brief, the Avenger PLS322-S-2 
model was determined to have a lower noise figure than the Star-Corn SR-3602 Mini, 
but the presence of electronics associated with a phase-locked local oscillator caused 
spurs in the output spectrum. Our software was unable to correct for rapid spurs, so we 
chose to continue using the Star-Corn SR-3602 Mini due to its history of stability and 
robust downconversion. The typical Y-factor between the beam at 8 degrees elevation 
and an absorber was 4.9 dB when using the Star-Corn SR-3602 Mini. Approximate 
estimates of noise contributions are 25 K from spillover, 50 K from the LNBF, 30 K 
from the atmosphere, and 5 K sky noise [4]. 

Radiation spillover efficiency was improved to the maximum theoretical value of 
91.6% by positioning the LNBF at an angle of 42° towards the dish (or 19.6 cm up 
from the bottom of the dish; see Appendix §8 for a more detailed calculation). We note 
that the efficiency can be boosted to 96% if the dish were 5.08 cm longer along the 
circumference, which could be achieved by machining a metal ring and sealing it along 
the circumference of the dish. We developed an LNBF shield (Figure 4) in order to 
further increase the spillover efficiency and lower the clamped back lobe of the LNBF. 
The shield was constructed from a Norpro 607 plastic casing funnel covered with tinfoil, 
and tested in the presence of the Haystack radorne. It was found to significantly reduce 
the radar interference, and increase the Y-factor by about 0.1 dB. Note, however, that 
we do not recommend using the shield for long-term spectrometer deployment, as it may 
fill with snow and ice in the winter or attract insects and birds. Alternative methods to 


2. Equivalently -0.3125 and +0.3125 in the I/Q spectrum. 




Figure 4: Antenna with LNBF shield. 


protect the spectrometer from radio-frequency interference include: (a) mounting the 
dish on a wall that prevents the LNBF from being in the line of sight of the Clarke 
belt 3 , (b) extending the size of the dish with ad-hoc metal petals, or (c) inverting the 
dish and keeping the LNBF shield so that water or snow cannot collect around the 
tinfoil. Lastly, we note that further software improvements include a high processing 
efficiency (70% of the theoretical maximum) with minimum time lost to calibration and 
calculations. Additionally, the power consumption of the two-channel unit was reduced 
from 60 W to 15 W due to switching our computer system from a Linux tower to the 
Intel NUC. 


3 Single-Channel Velocity Data 

3.1 Instrumentation 

To acquire velocity data, we used a network (Figure 5) of three single-channel ozone 
spectrometers located at the MIT Haystack Observatory in Westford, Massachusetts 
(42.618° N, 71.498° W), Union College in Schenectady, New York (42.817° N, 73.928° 

W), and Chelmsford High School in Chelmsford, Massachusetts (42.619° N, 71.367° 

W). The spectrometers are pointed at azimuths 190°, 180°, and 172°, respectively. 

Note that the beam physically intersects the ozone centered at 95 km above (37.953° 

N, 71.498° W), which is in the Atlantic Ocean about 200 km off the coast of Maryland. 

Data at Haystack were taken from January 24, 2009 through August 13, 2014. Data at 
Union College were taken from February 19, 2009 through June 8, 2015. Finally, data at 
Chelmsford High School were taken from December 25, 2011 through June 8, 2015. The 

3. The orbit at 35,786 km that is occupied by geostationary satellites and is a source of radio interference 
in the X and K u bands. 
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Figure 5: The geographic network of single-channel ozone spectrometers. The geotags, 
from left to right, indicate: Union College, MIT Haystack Observatory, and Chelmsford 
High School. 
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Figure 6: GPS calibration block diagram. 
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Figure 7: The physical implementation of the two-channel ozone spectrometer. 
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Figure 9: Nighttime trend in the velocity of ozone at 95 km. 
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six-channel and two-channel spectrometer readings were not included as we have yet to 
obtain one full year of data. Before processing the ozone data, we recalibrated the signal 
reference to the Haystack spectrometer, since its frequency calibration pulse generator 
circuit uses an atomic clock instead of a crystal oscillator. The atomic clock causes 
an improper calibration in the local oscillator of USB TV dongles. Using the setup 
in Figure 6, we first turned the calibration program on, and then initialized the GPS 
receiver for a stable 10 MHz reference signal, which was used to drive to the synthesizer 
to produce an output of 11.072 4545 GHz. We verified that the crystal oscillator has 
the same effect. 

3.2 Results and Discussion 

Figure 8 and Figure 9 show velocity data from the spectrometers at Haystack, Union, 
and Chelmsford. By convention, negative velocities denote winds that travel from north 
to south (northerly wind), and positive velocities denote winds that travel from south to 
north (southerly wind). For the general seasonal trend (Figure 8), data were obtained 
by binning over nighttime hours for each day of each year and plotting the weighted 
arithmetic mean for each month. The weighted standard deviation of the errors in each 
bin were used to plot the errorbars on each data point. Specifically, for each bin of size 
N containing velocity data points Vi and error data points Ei, the weights are defined 
to be 

Wi = ~2 . (5) 

The weighted arithmetic mean of the velocities is given by 

N 

v* = ^2 Vj - Wi ( 6 ) 

i=1 

and the error in each bin is taken to be the weighted standard deviation, which is the 
unbiased weighted estimate of the covariance matrix: 


N 

J2 w i\vi-v*\ 2 . (7) 

2=1 


We note that the meridional winds are northerly in the summer and southerly in the 
winter. The velocity appears to cross the zero mark twice, roughly around April and 
September, with northerly winds in between, and southerly winds otherwise. The max¬ 
imum amplitude in the average trend is 30 ms -1 , considering a minimum velocity of 
about -20 ms -1 in June-July and a maximum velocity of about +10 ms -1 in January. 

A nighttime trend is observed (Figure 9) in which there are southerly winds from -5 
through +1 hours local solar time, and northerly winds from +2 through +5 hours. The 
transition point appears to be within two hours after midnight, and the full transition 
is from about -20 ms -1 to +20 ms -1 in amplitude. Daytime data is not available due 
to the annihilation of ozone by ultraviolet photons from the Sun. Data were binned 
by day of the year, for each year, and averaged for each hour. The weighted standard 
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Figure 10: Comparison of the average meridional velocities of ozone and the oxygen green 
line by season and local solar time. 


deviation in the errors in each bin were again used to plot the errorbars for each data 
point. The data and trends from all three sites are in good agreement. Additionally, 
we note that magnitudes in the nighttime trend appear to be slightly larger than the 
magnitudes in the seasonal trend. 


4 Comparison with Millstone Hill Fabry—Perot Data 

4.1 Instrumentation 

The Millstone Hill High-Resolution Fabry-Perot Interferometer (FPI) is located near 
the MIT Haystack Observatory, specifically at (42.62° N, 71.45° W) where the mean 
local solar time lags behind universal time (UT) by 4 hours and 46 minutes. The 
instrument has a 100 mm aperture, and is pressure-tuned. A standard observation at 
30° elevation involves five optical measurements at azimuths 0°, 45°, 135°, 225°, and 
315° [15]. One of the measured parameters is the Doppler shift of the 557.7 nm “green” 
spectral line nightglow from atomic oxygen, which is a constituent of the aurora borealis 
[16] and has peak emission at an altitude of 95 km with a half-width half-intensity of 
8 km [17] [18]. This makes the altitude of the green line consistent with the altitude of 
the ozone line at 95 km. 

4.2 Dynamics of the Green Line 

By Equation 4, the Doppler shift of the green line produces information about the 
velocity of atomic oxygen in the upper mesosphere, which can be physically interpreted 
using winds and tides. The spectral line is the physical result of the 0( 1 S- 1 D) quantum 
energy level transition in the oxygen atom, and the Barth Mechanism plays the major 
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role in the nighttime green line emission [17] [19]. It is a two-step process: 


0( 3 P) + 0( 3 P) + M —)• O 2 + M, 

o( 3 p) + o;^o( 1 s) + o 2 , 


( 8 ) 

(9) 


first involving a three-body collision between two oxygen atoms and an atmospheric 
chemical M to produce excited diatomic oxygen, and then a two-body collision to pro¬ 
duce the excited state of atomic oxygen. 

4.3 Results 

Due to high optical background illumination from the solar excitation of atomic oxygen 
[20], the FPI does not take measurements during daytime hours. Scattering from other 
atmospheric particles and micro-meteors also hampers the ability to acquire precise 
measurements during the day. Therefore, in our analysis of the FPI data, we only con¬ 
sider the atomic oxygen nightglow between ±5 hours local solar time. The FPI wind 
data came from 692 files downloaded from the CEDAR 4 Archival Madrigal Database. 
The time period was chosen to be 1 January 2012 through 31 December 2015 (automat¬ 
ically pointing to the current day of the year, since it is not yet December). A Python 
script extracted and processed the data to produce the plots in Figure 10. 

Variations in the green line meridional velocity by season and local solar time appear 
to correlate well with respective trends in the ozone line. It is observed that the green 
line winds originate from the north in the summer of the northern hemisphere and from 
the south in the winter. Furthermore, there is a gradual increase from northerly to 
southerly winds between -4 and +5 hours local solar time. The uncertainties given by 
the FPI data are entirely statistical in nature, and are not associated with any flaws in 
the interferometer itself or the experimental procedure [15]. As a whole, the individual 
error values are low (less than 1.0 ms -1 on average), resulting in smaller errorbars 
compared to ozone. 

4.4 Consensus 

The ozone and green line velocity measurements come from two separate instruments, 
and the fact that they seem to correlate well indicates that these small-scale observations 
of the middle atmosphere may have some information about the dynamics of global 
winds. In particular, our results predict transitions in the direction of meridional winds 
by season and local solar time. 


4. Coupling, Energetics, and Dynamics of Atmospheric Regions. The database is located at http:// 
madrigal.haystack.mit. edu/. 
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5 Comparison with Other Data and Models 

5.1 Ascension Island Meteor Radar and Arecibo ISR 

We have also compared our data with averaged mesospheric horizontal winds measured 
over Ascension Island (8° S, 14° W) using very high frequency (VHF) meteor radar 
(MR). The measurements were taken from January 2002 through June 2011 at 95 km 
[8]. The results are presented as a color plot (Figure 12) and extracted by season (Figure 
11 ). Although the study was undertaken in a tropical latitude, the data clearly show 
a seasonal correlation in the amplitude of the meridional velocity between the meteor 
radar, FPI, and ozone data. Additionally, although the following is outside the altitude 
range of the MOSAIC spectrometer, a similar trend in meridional velocity is seen from 
Incoherent Scatter Radar (ISR) at the Arecibo Observatory in Puerto Rico targeting 
the F-region of the ionosphere at 150 km [21]. 

5.2 2D-SOCRATES, CIRA-72, HWM, and PWM 

The 2D-SOCRATES model is a theoretical model whose wind velocities are the result of 
a numerical solution to a partial differential equation with middle-atmosphere boundary 
conditions. The model is populated with a mix of chemical species, mostly hydrocarbons 
and chlorofluorocarbons [22] [23] . Therefore, there are no oxygen species present in this 
model to compare with the ozone and green line data. However, the meridional velocities 
at 95 km (Figure 13) have a similar trend in season, but with smaller magnitudes. 
This may be the result of zonally-averaged residual meridional wind calculations, or 
atmospheric tide extraction over all latitudes and longitudes, which would decrease the 
overall magnitudes. A similar trend in meridional velocity by season is also seen in the 
CIRA-72 wind model (Figure 14) near 90 km altitude at 43° N altitude [24], 

The Horizontal Wind Model (HWM) is an empirical model based on zonal and 
meridional wind data from numerous studies of the atmosphere under a range of clima- 
tologically averaged geophysical conditions [25] [26]. The original 1990 model, HWM90, 
was revised three times to produce HWM93, HWM07, and HWM15 in 1993, 2007, 
and 2015, respectively. While experimental data from certain regions of the atmo¬ 
sphere (notably the troposphere) seem to match well with the model, the team that 
produced HWM93 directly indicated significant discrepancies regarding wind data near 
the mesopause [25] . Other studies have cited inconsistencies between HWM and exper¬ 
imental observations, notably FPI measurements [9] [27]. 

Finally, we considered the 2D prevailing wind model, which incorporates climatolog- 
ically averaged data from 46 ground-based medium-frequency (MF) radar and meteor 
radar (MR) stations, as well as the HRDI satellite [28] . The seasonal trend in the merid¬ 
ional wind velocity in prevailing wind model did not correlate with data from ozone, 
Millstone Hill, or Ascension Island. 
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Figure 11: A comparison of the meridional wind velocities of the ozone, oxygen green line, 
and atmospheric winds from meteor radar. 


18 


































J FMAMJ JASOND 

Figure 12: Meridional mean wind from meteor radar measurements over Ascension Island. 
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Figure 13: Meridional wind velocities as given by the SOCRATES model. 
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Figure 14: Meridional wind velocities (black dots) as given by the CIRA-72 model. 
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6 Preliminary Observations of Ozone Temperature 

In previous work, the theoretical ozone temperature in the mesosphere and lower ther¬ 
mosphere (MLT) was determined to be around 170 K due to particle collisions and 
ultraviolet photon absorption. Measurements of the ozone line by initial Very Small 
Radio Telescope (VSRT) hardware had insufficient spectral resolution to provide an ac¬ 
curate measurement of the ozone temperature, but did suggest that the temperature was 
below 180 K [29]. Recent nighttime data taken from the network of three ozone spec¬ 
trometers (Haystack, Chelmsford, Union) from 20 January 2015 through 10 May 2015 
(winter and spring) reveal an average ozone temperature of 170±12 K. However, looking 
at one individual measurement from only the Haystack spectrometer reveals an ozone 
temperature of 120 ±118 K. Thus, singular measurements of the ozone temperature 
seem to have erratically high errorbars, while the average has a lower overall tempera¬ 
ture. By comparison, satellite observations of the oxygen green line show temperatures 
in the 180-200 K range at 95 km altitude and 45° N latitude [17]. Further observations 
by Rayleigh Mie-Raman lidar and potassium resonance lidar show a seasonal temper¬ 
ature variation of 160-190 K at 95 km [6]. Lastly, sodium lidar measurements in Fort 
Collins, Colorado (41° N, 105° W) show a seasonal temperature variation of 180-190 K 
at 95 km [30]. 

7 Future Work 

7.1 The Global MOSAIC Network 

In contrast to satellite-based instruments, ground-based spectrometers have the advan¬ 
tage of time resolution but are weak in spatial resolution [28]. In order to achieve a 
greater spatial coverage in latitude and longitude, M. J. Kosch of Lancaster Univer¬ 
sity has proposed to expand the MOSAIC network to encompass eighteen two-channel 
spectrometers that will be located across the United States of America (Table 1), Eu¬ 
rope (Table 2), Africa (Table 3), and Antarctica (Table 4). Potentially, some of these 
spectrometers would also be pointed toward the east or west in order to obtain data 
for the zonal wind velocity of ozone at 95 km. This would provide a complete map of 
the two-dimensional motion of the ozone, and also reveal more about the global wind 
dynamics at that altitude. Increased averaging in the data (by increasing the number 
of instruments) would also fine-tune the resolution of ozone temperature. 

7.2 Further hardware upgrades 

Ideally, the Global MOSAIC Network would make use of the atomic clock-based fre¬ 
quency calibration pulse generator, instead of the one that is based on the crystal 
oscillator. The MIT Haystack Observatory has placed an order for several new atomic 
clocks that will eventually be sent to M. J. Kosch for field deployment. 

7.3 Helping to resolve the theory of pole—to—pole flow 

Lastly, we note that there exists a theory of inter -hemispheric wind circulation in the 
mesosphere, also known as pole-to-pole flow [12]. As shown in Figure 15, there is a 
seasonal trend in the mesospheric wind, in which air from the south pole moves upward 
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Table 1: The American MOSAIC Network 


Spectrometer Location 

Town 

State 

Latitude (°N) 

Longitude (°E) 

Chelmsford High School 

Chelmsford 

Massachusetts 

42.62 

288.63 

MIT Haystack Observatory 

Westford 

Massachusetts 

42.50 

288.53 

Bridgewater State University 

Bridgewater 

Massachusetts 

41.99 

288.99 

Lynnfield High School 

Lynnfield 

Massachusetts 

42.54 

288.97 

Union College 

Schnectady 

New York 

42.80 

286.07 

University of North Carolina 

Greensboro 

North Carolina 

36.07 

280.17 

Alaska Pacific University 

Anchorage 

Alaska 

61.19 

210.20 


Table 2: The European MOSAIC Network 

Spectrometer Location 

Region 

Country 

Latitude (°N) 

Longitude (°E) 

Ny-Alesund 

Svalbard 

Norway 

78.93 

11.92 

Sodankyla 

Lapland 

Finland 

67.42 

26.58 

University of Lancaster 

Lancashire 

United Kingdom 

54.05 

357.20 


Table 3: The African MOSAIC Network 

Spectrometer Location 

Region 

Country 

Latitude (°N) 

Longitude (°E) 

Cairo 

Cairo Governorate 

Egypt 

30.05 

31.23 

Bahir Dar 

Amhara 

Ethiopia 

11.60 

37.38 

Nairobi 

Nairobi County 

Kenya 

358.72 

36.82 

Windhoek 

Khomas 

Namibia 

337.43 

17.08 

Sutherland 

Northern Cape 

South Africa 

327.61 

20.66 


Table 4: The Antarctic MOSAIC Network 


Spectrometer Location 

Region 

Continent 

Latitude (°N) 

Long (°E) 

Halley Bay (UK) 

Brunt Ice Shelf 

Antarctica 

284.42 

333.43 

SANAE IV (South Africa) 

Vesleskarvet 

Antarctica 

288.33 

357.16 

Amundsen-Scott Station (USA) 

South Pole 

Antarctica 

270.00 

0.00 
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Figure 15: Inter-Hemispheric Circulation. 


through the lower thermosphere in the summer, and remains in the lower thermosphere 
until the winter, when it moves back into the mesosphere. This may be related in some 
form to the seasonal change in the polarity of the meridional winds as observed by 
the single-channel ozone spectrometers, FPI, and meteor radar. Earlier work by Yuan 
indicated a summer-to-winter pole-to-pole meridional flow of 17 ms -1 at 86 km [30], 
which is within the range of magnitudes that we have observed at 95 km. 

Using the Global MOSAIC Network, further studies of the meridional and zonal veloc¬ 
ities of ozone could reveal novel information about the Earth’s pole-to-pole circulation 
as well as the dynamics of air flow between the MLT and other regions of Earth’s atmo¬ 
sphere. In conclusion, we hope to see the two-channel hardware replicated and deployed 
at key locations around the world in the coming years in order to realize a global net¬ 
work of spectrometers that will provide higher-resolution data about the dynamics of 
ozone in the middle atmosphere. 


8 Appendix: Optimizing Antenna Geometry 

8.1 Summary 

Here, we develop a theoretical model for the antenna geometry, which is incorporated 
into the physical implementation of the two-channel ozone spectrometer. A set of 
equations are developed that relate the angle of the low-noise block downconverter 
feedhorn antenna (LNBF) to the radiation flux that hits the parabolic reflector. A 
simulation in Python determines that the spillover efficiency is maximized, and hence 
power loss is minimized, when the LNBF is pointed towards the dish at an angle of 41 
degrees. 
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8.2 
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Figure 16: Parabolic antenna geometry. 

In Figure 16, we model the Winegard DS-4040 dish as an offset parabola with 0.4572 
m diameter that extends into a region spanned by the {x, y, z} basis with origin 0\ and 
focus / 2 = 10.0 cm. The dish itself is the solid black line in the positive octant, while 
the full parabola includes the dashed line extending into the negative octant. The unit 
vector x represents East, y represents North, and z represents Up. In unit spherical 
coordinates, the location of an arbitrary point in this basis is given by 

{ x = cos 9 cos <p 

y — cos 9 sin </2 (10) 

z = sin 9 


where 9 and ip are the altitudinal and azimuthal angles, respectively. The antenna feed 
is located at origin 0 2 with focus / = 10.6 cm 5 , and any vector emerging from this 
origin is described by the “primed” basis, which consists of the unit vectors {x',y',z'}. 

As in Figure 16, to rotate the antenna by some angle 7 , we choose to keep the x'-axis 
constant and then apply a rotation matrix on the j/z-plane: 

w 

z' 

This results in the following transformed coordinate values: 

! x' = x 

y' = ycosy + zsiny ( 12 ) 

z' = z cos 7 — y sin 7 . 

5. The beam focus center was measured to be 0.6 cm into the circular waveguide of the LNBF, and not 
entirely at the “dot” on the cover of the feed. 


COS 7 sm 7 
- sin 7 cos 7 


( 11 ) 
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The transformed LNB angle 7 ' is: 


7 / = arctan 


(x' 2 + z’ 2 ) 


1/2 


(13) 


The location of some point p is described by the parabola equation 

p 2 =x 2 + z 2 = 4/a, (14) 

where a is the projection of p onto the negative y-axis. However, since we want to look 
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Figure 17: Coordinate scaling. 


at the situation from the perspective of the antenna, we must shift over to origin O 2 . 
As depicted in Figure 17, we take the unit coordinates emerging from 0 2 and choose 
some scaling factor d that makes p = (x'd) x' + (y'd) y' + ( z'd ) z'. Consequently, 

h - y’d = a. ( 15 ) 

and the parabola equation is modified accordingly: 

(x /2 + y ,2 )d 2 = 4f(f 2 -y'd). (16) 


We realize that this is a quadratic equation in d with coefficients 


a = (x' 2 + z' 2 )/ (4/) 
b = y 

c= -h 


and solution 


d = 


—6 + sjd 2 — 4ac 
2a 


(17) 


(18) 


where we only consider the positive square root as it is physically realizable. The offset 
parabola, in the scaled coordinate system, is now extended into circle with radius r and 
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a j/-offset of r to account for the antenna actually being shifted by that amount. The 
quantity R is defined to be 


R = 



r) 2 + ( x'd ) 2 

j.2 


(19) 


and is used as follows: if R = 1 , then the beam emanating from the antenna is at the 
edge of the dish. If R < 1, then the beam spills over past the edge (i.e. it is not on 
the circle), resulting in a power loss because not all of the radiation is striking the dish. 
Quantitatively, the spillover efficiency e s is the ratio of flux intensities [31] given by 


JJ dfl G(9 , <p, 7 ) cos 6 


{r<l} 


dfl G(9,ip, / y) cos9 


( 20 ) 


where f l is the solid angle subtending the dish (for all r < 1 in the numerator, and 
across the entire region in the denominator) and G(9,ip, 7 ) is the gain function of the 
antenna, which we define as a Gaussian: 

G{0, ip, 7 ) = exp (-feV 2 ) ( 21 ) 

= exp (—fe arctan /) ( 22 ) 

where 

cos 7 cos 9 sm ip + sin 7 sin 0 

The quantity arctan / is the antenna angle in the O 2 basis and fe = —0.692/(22) 2 is a 
normalization constant that takes into account the half-width half-maximum (HWHM) 
power at 22 degrees orientation. 


8.3 Results 

Due to a lack of efficient computational resources, direct symbolic integration of Equa¬ 
tion 20 is not feasible with Wolfram Mathematica, so we perform numerical integration 
in Python. A plot of the output of the program is shown in Figure 18. We find that 
a maximum efficiency of 0.916 occurs when the antenna is pointed at 7 = 41 degrees 
from the horizontal. 


9 Appendix: LNBF Characterization 

9.1 Introduction 

The Y-factor is widely used to measure the noise temperature of an amplifier. It is 
defined to be the ratio of the amplifier power when “hot” to the power when “cold,” 
and mathematically, this is represented by the following equation: 

y. _ Phot _ Tirnb T TlNBF (241 

Pcold Tsky + TlNBF 
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Figure 18: Theoretical spillover efficiency versus antenna angle. 


where the “hot” power, 


-Phot = T am b + Ilnbf, 


(25) 


is the power measured when the LNBF is pointed at an ambient source with temperature 
T am b ~ 300 K, and the “cold” power, 


-P co id = T sky + Ilnbf, ( 26 ) 

is the power measured when the LNBF is pointed at the sky with temperature T sky ~ 
3 K + k. Here, the 3 K is due to the background radiation in the sky, and k is some 
constant (in units of Kelvin) that characterizes the attenuation of the microwave signal 
due to atmospheric absorption. 

9.2 Experimental Procedure 

We use a C program to continuously measure and output the power values of the 
amplifier in dB. Using the results from Appendix §8, we point the LNBF at 42 degrees 
toward the parabolic reflector. The dish itself is then pointed 6 42 degrees toward the 
sky. This produces a measurement of P co id in the program output. Next, a Cuming 
Microwave pyramid-shaped carbon fiber microwave absorber is wrapped around the 
feedhorn of the LNBF. Specifically, one pyramid is cut and placed into a styrofoam cup, 
which can be easily fastened over the LNBF (Figure 19). The program then outputs 
a measurement of Phot- Since the power units are in log-scale, the Y-factor can be 
obtained by taking the absolute value of the difference of the two measurements: 

1 [dB] = abs(dB a hsorb er dB no absorber )■ (27) 
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Status 

LNB Model 

Channel 1 (dB) 

Channel 2 (dB) 

Without absorber 

Avenger 

PLL322-S-2 

-46.5 

-41.0 

With absorber 

-40.1 

-35.2 

Y-factor 

6.4 

5.8 

Without absorber 


-40.8 

-36.5 

With absorber 

Chaparral 

-35.6 

-33.2 

Y-factor 


5.2 

3.3 

Without absorber 


-42.5 

-39.1 

With absorber 

Octagon 0TL50 

-38.4 

-34.4 

Y-factor 


4.1 

4.7 

Without absorber 

Star-Corn 

-39.9 

-34.2 

With absorber 

SR-3602 Mini 

-37.1 

-30.7 

Y-factor 

2.8 

3.5 

Without absorber 

Avenger 

KSC322-2 

-39.2 

-37.6 

With absorber 

-34.4 

-33.1 

Y-factor 

4.8 

4.5 


Table 5: LNBF Y-Factors 


LNB Model 

Channel 1 (K) 

Channel 2 (K) 

Avenger PLL322-S-2 

37.8 

44.4 

Chaparral 

52.9 

109 

Octagon 0TL50 

76.9 

62.0 

Star-Com SR 3602 Mini 

143 

99.0 

Avenger KSC322-2 

60.6 

66.4 


Table 6: LNBF Noise Temperatures 


LNB Model 

Channel 1 (dB) 

Channel 2 (dB) 

Avenger PLL322-S-2 

0.53 

0.62 

Chaparral 

0.73 

1.4 

Octagon 0TL50 

1.0 

0.84 

Star-Com SR 3602 Mini 

1.7 

1.3 

Avenger KSC322-2 

0.82 

0.90 


Table 7: LNBF Noise Figures 
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Figure 19: The two-channel ozone spectrometer with an attached microwave absorber 
during Y-factor testing. 



SURFACE WATER-VAPOR DENSITY, p (gm/m 3 ) 


Figure 20: Atmospheric zenith opacity (right-hand vertical axis) and atmospheric decibel 
attenuation (left-hand) axis versus surface water-vapor density. Note that while Meeks 
presents data for a 15 GHz signal, the ozone spectrometer is actually tuned to the 11 GHz 
spectral line. 

9.3 Results 

Measurements and computed Y-factors for five different LNB models are given in Table 
5. We proceed by calculating the signal attenuation constant k. From [32], we find that 
the saturation pressure of water vapor is given by the equation 

exp (77.3450 + 0.0057T - 7235T" 1 ) , 

Pws = ^- 1 (28) 

where T is the dry bulb temperature of moist air. Using a mean temperature of 290.9 
K from weather records for June 26, 2015, we compute a saturation pressure of p ws = 
2026.79 Pa [33]. The water vapor density (in kg • m~ 3 ) is given by the equation 

p w = 0.0022 x ^ (29) 

6. In practical runs of the ozone spectrometer, the dish should actually be pointed to around 9 degrees 
toward the sky. The fact that the dish is pointed 42 degrees toward the sky in this experiment is attributed 
to human error. Note, however, that this does not affect the resulting Y-factor measurements. 


28 










Figure 21: Zenith correction geometry. 


and given T = 290.9 K, we compute a water vapor density of 15.33 g • nr -3 [32]. 
According to Figure 20, a water vapor density of 15.33 g • m -3 corresponds to an 
atmospheric zenith opacity ~ 0.03. The temperature of the atmosphere follows as 

^atmosphere = K = 0.03 X T am b ~ 8 K. (30) 

This yields T s k y = 3K + 8K = 11K. However, Figure 20 cites the zenith opacity, and 
the dish is not pointed at zenith. If the atmosphere has a height h and our dish is 
pointed at some angle 9 = 42° , then the length £ that forms the hypotenuse of the right 
triangle in the direction of the dish (see Figure 21), is given by 

£ = hcsc8. (31) 


We would expect a similar trend in the temperature: 


rpl 

atmosphere 


^atmosphere CSC (42 ) 


12 K. 


(32) 


Thus, in reality, Tsky = 3 K +12 K = 15 K. Now that we have determined k, we return 
to Equation 24 to compute the LNB noise temperatures. The results are presented in 
Table 6. Finally, we convert the noise temperature to the noise figure (NF) in decibels, 
which is defined as 


NF [dB] 


101og 10 


/ 290 K + r LNB \ 
v 290 K ) ' 


(33) 


The noise figures for all six LNB models are presented in Table 7. 


9.4 Consensus 

The theoretical noise temperatures for the Star-Corn SR-3602 Mini are consistent with 
those measured for the six-channel ozone spectrometer, which was developed by A. E. 
E. Rogers and E. True [3]. This agreement ensures that our calculations are robust. 
The LNB with the best noise performance is the Avenger PLL322-S-2, which has a 
noise figure of 0.53 dB on Channel 1 and 0.62 dB on Channel 2. In contrast, the LNB 
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with the worst noise figure appears to be the Star-Com SR 3602 Mini , with 1.7 dB on 
Channel 1 and 1.3 dB on Channel 2. However, upon testing the Avenger PLL322-S-2 
in the field, we discovered that there was a significantly higher frequency drift in the 
Doppler shift of the ozone line at 11 GHz. This is due to a phase-locked local oscillator 
inside the Avenger LNB circuitry, which makes it much less stable for long-term data 
acquisition as our C program is unable to calibrate the LNB frequency at faster rates. 
In conclusion, we have chosen to continue using the Star-Com SR-3602 Mini as our 
LNB of choice, largely due to its past history of being a stable antenna feed for the 
ozone spectrometer hardware. 


10 Appendix: Code 

Most of the source files used in this project are publicly available at the ozone spec¬ 
trometer GitHub repository, which is located at bit. ly/lOsnYf 8. The Python script 
used to generate the file download from CEDAR Madrigal (§4), as well as the file log, 
is available here: bit. ly/lVL21wS. 
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